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production mechanism dominates over the fragmentation production mechanism at the all relevant 
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I. INTRODUCTION 



The phenomenology and the theory of processes involving the production of heavy 
quarkonia have been vigorously developed for last decade following the measurement of 
the transverse- momentum (py) spectra of prompt J/ip mesons by the CDF Collaboration 
at the Tevatron Collider [l| and by the ZEUS and HI Collaboration at the HERA Collider 
{2I, 3)]. The color-singlet model previosly proposed to describe a non-perturbative tran- 
sition of a QQ pair to final-state quarkonium was extended in a natural way within the 
formalism of nonrelativistic QCD (NRQCD) [5j. The extended model takes into account 
QQ-pair production not only in the color-singlet state but also in the color-octet state. 

The NRQCD formalism makes it possible to calculate consistently, by perturbation theory 
in two small parameters (strong coupling constant a s at the scale of the heavy-quark mass 
and the relative velocity v of quarks in quarkonium), not only the parton cross sections 
for quarkonium production processes through the fusion of Q and Q quarks but also the 
universal fragmentation functions for parton splitting into various quarkonium states 

Usually it is assumed that, in the region of large quarkonium transverse momenta (p? ^> 
M, where M is the quarkonium mass), the fragmentation production mechanism is more 
adequate, because it takes into account effectively high order corrections in a s , than the 
mechanism associated with the fusion of a heavy quark and a heavy antiquark produced in 
a hard subprocess, which is calculated in the lowest order (LO) in strong coupling constant 
a s . However, the numerical value of the critical transverse momentum is unknown a'priori. 

In the energy region of the Tevatron and HERA Colliders the main contribution to the 
heavy quarkonium production cross sections comes from the gluon-gluon or the photon- 
gluon fusion at the small values of the argument x of the gluon distribution function. In 
the conventional collinear parton model [71, the initial-state gluon dynamics is controlled 
by the Do^-G.bov-Lipatov-Alta.em-P.isi (DGLAP) e TOlu «on e qU a tl on g la tHs 
approach, it is assumed that S > p 2 ^> Aq CD , where yS is the invariant collision energy, 
/i is the typical energy scale of the hard interaction, and Aqcd is the asymptotic scale 
parameter of QCD. In this way, the DGLAP evolution equation takes into account only one 
big logarithm, namely / Aqcd) ■ In fact, the collinear-parton approximation is used, and 
the transverse momenta of the incoming gluons are neglected. 

In the high-energy limit, the contribution from the partonic subprocesses involving t- 
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channel gluon exchanges to the total cross section can become dominant. The summation 
of the large logarithms ln(v#//f) in the evolution equation can then be more important 
than the one of the ln(/// Aqcd) terms. In this case, the non-collinear gluon dynamics 
is described by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation [9|. In the 
region under consideration, the transverse momenta (fcr) of the incoming gluons and their off- 
shell properties can no longer be neglected, and we deal with reggeized t-channel gluons. The 
theoretica 
approach 



frameworks for this kind of high-energy phenomenology are the /c^-factorization 
10, 11] and the (quasi)-multi-Regge kinematics (QMRK) approach 3]. Last 



one is based on effective quantum field theor y im plemented with the non-abelian gauge- 
invariant action, as suggested a few years ago [141 ] . However, the /c^-factorization approach 
has well-known principal difficulties 15|] at next-to-leading order (NLO). By contrast, the 
QMRK approach offers a conceptual solution of the NLO problems 161 ] . 



Our previous analysis of J / ip meson production at the Fermilab Tevatron Colliders using 
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191 ] has shown the 



the high-energy factorization scheme and the NRQCD formalism fl7. 
dominant role of the color-octet intermediate state 3 S[ 8 ^ contribution at the large px re- 
gion in the fusion production and in the fragmentation production with the gluon splitting. 
We found also the good agreement between the values of non-perturbative matrix element 
(NME) (O j ^[ 3 S{ 8) ]) obtained by the fit in collinear parton model and in the QMRK ap- 
proach. The similar analysis for the J/ip photoproduction at the high energy is the subject 
of study in this paper. The situation should be different because of the c— quark fragmen- 
tation via the color-singlet state dominates over the gluon fragmentation via the color-octet 
state in J/tp photoproduction. 

This paper is organized as follows. In Sec. [Til we discuss basic ideas and formulas of the 
fragmentation approach, including the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations for the c— quark and gluon fragmentation functions into J/ ip meson. In 
Sec. IIHI we present relevant LO squared amplitude in the QMRK approach and master 
formulas for the differential cross sections. In Sec. IIVI we discuss our results for the J/ip 
photoproduction at the HERA Collider. In Sec. [V] we summarize our conclusions. 
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II. FRAGMENTATION MODEL 



Within the framework of the NRQCD, the fragmentation function for quarkonium H 
production can be expressed as a sum of terms, which are factorized into a short- distance 
coefficient and a long distant matrix element Q]: 

D(a^H) = Y J D(a^QQ[n})(0 H [n}). (1) 

n 

Here the n denotes the set of color and angular momentum numbers of the QQ pair, in which 
the fragmentation function is D(a — > QQ[n\). The last one can be calculated perturbatively 
in the strong coupling constant a s . The non-perturbative transition of the QQ pair into 
H is described by the NMEs (O n [n]), which can be extracted from experimental data. To 
LO in v, we need to include the cc Fock states n = 3 S i i\ 3 Sf\ ^^^Pf if H = J/ip(ip'). 
The recipes of calculations for heavy quarkonium fragmentation functions at the initial scale 
1^ — l^o — 2 m c i n LO NRQCD formalism are well known [6] and detailed formulas for the 
gluon and c— quark fragmentation via the different intermediate states are presented, for 
example, in Ref. 201 ] . 

To calculate the J /if) production spectra, we need the fragmentation functions at the 
factorization scale /i ^> M. Then, large logarithms in fi/M appear, which have to be 
ressumed. It is achieved by using the DGLAP equations 

j^D^ H (g,fi) = [ —Pba{-,ot s {n))D h ^ H {x,n)i (2) 

where P& a ( — , ct s (u)) are the timelike splitting functions of parton a into parton b. The initial 

X □ 
conditions D a -+H\£i Ho) f° r a — c ->9 an d H = J/if)(if)') are taken from Ref. [20(. The initial 

light-quark fragmentation functions are set equal to zero. They are generated at larger scale 

jj,. However, in the case of J /if) production at the HERA Collider their contribution is very 

small. Here, a = g,c, c. To illustrate the £ and /x dependence of the relevant fragmentation 

functions we show in Figs. [1] and [2] the c— quark and the gluon fragmentation functions at 

the different values of /j. We use LO in a s approximation for the splitting functions. The 

dashed curves are obtained neglecting the c — g coupling during the evolution (a = c or 

a = g), the solid curves are obtained using the system (j2J), which takes into account cross 

transitions between gluons and c— quarks during the evolution. 
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III. BASIC FORMULAS 



Because we take into account the gluon and c— quark fragmentation into the J/i[) mesons, 
we need to sum contributions from LO partonic subprocesses in which these ones are pro- 
duced. 

In the collinear parton model we take into consideration the following subprocesses: 

1 + 9 -> c + c, (3) 
7 + 9(9) 9 + Q(Q), ( 4 ) 

where q = u,d, s. The squared amplitudes of the subprocesses ([3]) and PJ are well known. 
In the case of QMRK approach we have only one dominant LO partonic subprocess: 

7 + R -> c + c, (5) 

which is described by the Feynman diagrams shown in Fig. [31 Here R is the reggeized gluon. 
Accordingly amplitude of the subprocess (j3J) can be presented as follows: 

A(j + R^c + c) = -ig s e c eT a U{p u m c )[ 7 Q . Pl - 9 ! 2 + ^ 2 7^ + 

(Pi - qi) - ™ c 

p Pi - + m c mc )^( gi ) (n -)^ ( 6 ) 

(Pi - 9ar - «i 2 

where = P^/E N , P£ = E N (1, 0, 0, 1) is the four-momentum of the proton, £jv is the 

energy of the proton, e a (q\) is the polarization four- vector of the photon. For any four- 
momentum we define k~ = (k ■ n~). Note, that reggeized gluon four-momentum can be 
written as 

(£ = x 2 E N {n-y + (£ r , (7) 

where q 2T = (0, q 2T ,0), q\ = q 2 T = —q£ T and x 2 is the fraction of the proton longitudinal 
momentum passed on to the reggeized gluons. 
The squared amplitude obtained from (jSJ) reads: 

647T 2 e 2 c a K a , 
I 2 !! 2 



\A{i + R - c + c)| 2 = \ s [2m 2 fe - up' 2 f + c£ T ut((pT) 2 + (p 2 -) 2 )], (8) 



where the bar indicates average (summation) over initial-state (final-state) spins and colors, 
u — u — m 2 c , t — t — ml, u—[q x — p 2 ) 2 , t — {q% — pi) 2 , a = e 2 /Air, a s = g 2 J^, e c = +2/3. 
The obtained result ([S]) satisfies gauge invariant conditions, 

lim \A(i + R^> c + c)| 2 = 0, 



and agrees with the fc^-factorization formula from Ref. 25( if we take into account relation 



between amplitudes in the QMRK and the k?— factorization approach 

(X 2 E; A2 



\Akt{i + R — » c + c 



-\|2 



Q2T 



■\Ain + R 



c + c)\ 2 . 



(9) 



At HERA, the cross section of prompt J/ip production was measured in a wide range 
of the kinematic variables W 2 = (Pn + <?i) 2 , Q 2 = — Qi, Xi = (Pn ' Qi)/(Pn • k), z — 
(Pn • V)l (Pn ■ qi), and pr, where P N , q± = — /c' M , and p^ are the four-momenta 

of the incoming proton, incoming lepton, scattered lepton, virtual photon, and produced 
J/ip meson, respectively, both in photoproduction 3], at small values of Q 2 , and deep- 
inelastic scattering (DIS) 3J, at large values of Q 2 . At sufficiently large values of Q 2 , the 
virtual photon behaves like a point-like object, while, at low values of Q 2 , it can either act 
as a point-like object (direct photoproduction) or interact via its quark and gluon content 
(resolved photoproduction). Resolved photoproduction is only important at low values of z 
and pt- The subject of our study is large px photoproduction processes at the Q 2 ~ only. 

In the region of px 3> m c we consider c-quark and J / ip- meson as massless particles, so 
that the fragmentation parameter £ is related to the c— quark and J/ip— meson four-momenta 
as follows: 

P M = (10) 

We only used non-zero c— quark mass in the initial factorization scale /j,q = M — 2m c and, 
correspondingly, in the definition of the \i 



lp T + 4m 2 . 

Let us first present the relevant formula for the double differential cross sections of J/ip 
direct photoproduction in the collinear parton model via the partonic subprocess (J3j): 

da ep 



dzdpj- 



where 



x 2 



a f ( \ f 1 d ^n ic \\ M (l + 9 -> c + c)\ 2 x 2 G(x 2 ,n) 

^MJ^j^^) 1GKzi z- z){xiXtS) 2 < 



Pt 



x u 



Pt 



z + 



Pt 



ZXib(t, — Z) Zb{l — Z) X\Zb 

G(x 2 , n) is the gluon collinear distribution function in the proton, for which we use numerical 
CTEQ parameterization 21[ and / 7 / e (a;i) is the quasi- real photon flux. In the Weizacker- 
Williams approximation, the latter takes the form 



f-t/eiXl) 



a 
2^ 



1 + (1 



Xl 



O 2 
In max 



Xi 



Q 



2 
min 



2m 2 e Xi 



Qmin Q 



2 

max , 



(12) 
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where Q 



2 

min 



Q 



2 

max 



1 GeV 2 



m?xf/(l — x\) and Q max is determined by the experimental set-up, e.g. 

3- 

The formulae in the case of the partonic subprocess (jlj) looks like ( JTTT) . The J/ip mesons 
are produced also via decays of excited charmonium states and B— mesons. In fact, only the 
contribution of the decays ip' —>■ J/ipX is large at the relevant kinematic condition. In the 
approximation used here we can estimate this one simply: 



da(ep -> J/tpX) = da drrect {ep -> J/iffX) + da direct {ep -> ?//X)Br(^' J/^X), (13) 

where Br( , / — > J/ipX) is the inclusive decay branching fractions, da dtrect (ep — > can 
be obtained formally from da direct (ep — > J/ipX) after replacement ((9 J// ^[n]) — > (0^'[n]). 

In the QMRK approach the double differential cross sections of J/ ^ direct photoproduc- 
tion reads: 



da ep 
dzdp\ 



^— J d(j) 2 J J dxif l/e (xi) J yD c ^((,/j) x 



\(2ir)< J J q 2T 
{x 2 E N ) 2 $(x 2 ,c£ T ,ti) 



1-4(7 + R 



c + c 



(14) 



z(£ - 2)(xix 2 5') 2 

where 02 is the angle between q 2 T and the transverse momentum of the J ftp meson, 
$(x2, q^T) A 4 ) is t ne unintegrated reggeized gluon distribution function. The collinear and 
the unintegrated gluon distribution functions are formally related as: 



x 2 G(x 2 ,n) 



dql T $(x 2 ,ql T ,n), 



(15) 



so that the normalizations of Eqs. (fl4l) and (iTTj) agree. In our numerical calculations, we 
use the unintegrated gluon distribution function by Kimber, Martin, and Ryskin (KMR) 
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24) of the 



22( | . which gives the best description in the framework of QMRK approach 
heavy quarkonium production spectra measured at the Tevatron [l|, Q] ■ 

To compare predictions of the fragmentation and the fusion mechanisms we need to 
calculate J /if) photoproduction cross section using the fusion model too. In the case of 
the LO collinear parton model calculation we use squared amplitudes for the color-singlet 
contribution from Ref. Q and for the color-octet contribution from Ref. 26]. Working in the 



QMRK approach with the fusion model we use analytical results obtained in Refs. 
All numerical values of parameters are the same as in the Ref. 17]. 
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28] 
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IV. RESULTS 



We now present and discuss our numerical results. 

Let us compare J/ if) meson px spectra obtained using the fragmentation and fusion models 



whithin the QMRK approach. As it was shown recently [171 . |281 |. we can describe the 
experimental data from HERA Collider j^, jj] in the fusion model without the color-octet 
contribution, g.e. in the color-singlet model |4| (see the curve 1 in the Fig.H]). The agreement 
is good for the all measured J /if} transverse momenta up to p\ = 50 GeV 2 . The prediction 
of the fragmentation model is shown in Fig. H]for p\ > 10 GeV 2 only (the curve 2). We 
see that the curve 2 tends to the curve 1 as p\ increases but it is being lower at all relevant 
p\ till p\ = 10 3 GeV 2 . The result is the same as in the hadronic J/tp production 
although in the photoproduction the c— quark fragmentation via the color-singlet state 3 S^ 
is dominant, oppocite to the hadroproduction where the gluon fragmentation via the color- 
octet state 3 S[ 8 ^ is dominant. So, we see that in the QMRK approach the slope of the pt 
spectra sufficiently depends on the initial reggeized gluon transverce momentum, which is 
ordered by the unintegrated distribution function, and both production mechanismes, the 
fusion and the fragmentation, have similar p?— slopes at the large transverse momentum of 
the J/ if) meson. 

The situation in the relevant LO calculations within the collinear parton model was 
studied earlier in Refs. 
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271 ] . It was shown that at approximately p\ > 100 GeV 2 



the contribution of the c— quark fragmentation production exeeds the contribution of the 
fusion mechanism with the color-singlet LO amplitudes. Note, that in both of the discussed 
approaches, the fragmentation and the fusion within the CSM in LO, do not describe data. 
However, if we take into account the color-octet contribution, with the values of the color- 
octet NMEs fixed at the Tevatron, we see that the fusion model describes data well, as it is 
shown in Fig. [5] (the curve 1). The contribution of the fragmentation model increases too, 
because of the additional piece coming from the gluon fragmentation via 3 S[ 8 ^ state in the 
partonic subprocesses 7 + q — > g + q (the curve 6 in Fig. [5]). Moreover, at p\ > 500 GeV 2 
the gluon fragmentation dominates over the c— quark fragmentation (the curve 5), but their 
sum (the curve 4) is still less than the fusion model prediction in the NRQCD approach up 
to p\ = 2 ■ 10 3 GeV 2 . 

For completeness, we present here the z— spectra of J/if> calculated in the QMRK ap- 
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proach (Fig. [6]) and in the collinear parton model (Fig. [7]) with the experimental data from 
HERA Collider at the px > 3 GeV 3[ . We see that in both approaches the fusion mechanism 
is dominant for all z, excepting the region of very small z < 0.1. It is well known that at 
small z the resolved J/ip meson photoproduction becomes important. However, it is evident 
that in the resolved photoproduction the relation between the fusion and fragmentation 
mechanisms will be the same as in the direct photoproduction. 

V. CONCLUSIONS 

It was shown that working within NRQCD both in the collinear parton model and the 
QMRK approach, we have predicted the dominant role of the fusion mechanism as com- 
pared to the fragmentation one in the J/ip meson photoproduction at all relevant values 
of charmonium transverse momenta. Therefore, it is impossible to extract the c— quark or 
gluon fragmentation function into J/ip meson from the photoproduction data similar as in 



the hadroproduction processes, as it was shown previosly [18|]. Opposite to the assumption 
in the collinear factorization, that the fragmentation contribution is enhanced by powers of 
{jPr/wfy relative to the contribution of the fusion model, we have demonstrated that such 
an enhancement is absent in the non-collinear factorization approach, QMRK. We have ob- 
tained equal slopes of the px— spectra for both production mechanisms at the asymptotically 
large pt in the QMRK approach. 
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FIG. 1: The fragmentation function D c -*J/il>{ z i M 2 ) a ^ the fJ 2 = 10 GeV 2 (the curves 1 and 3) and 
/i 2 = 300 GeV 2 (the curves 2 and 4). The solid lines are results of Eqs. (|2j), the dashed are obtained 
neglecting the c — g coupling. 
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FIG. 2: The fragmentation function Dg^j/^^z, n 2 ) at the fi 2 = 10 GeV 2 (the curves 1 and 3) and 
fi 2 = 300 GeV 2 (the curves 2 and 4). The solid lines are results of Eqs. (|2|), the dashed are obtained 
neglecting the c — g coupling. 
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7, Qi Q, pi 7, Qi Q, pi 

v/\/\/\/> 

--^-^qqqqI 



FIG. 3: The Feynman diagrams of the process 7 + R — > c + c 
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FIG. 4: pt spectrum of prompt J/ip in ep scattering with E p = 820 GeV, E e = 27.5 GeV, 60 GeV 
< W < 240 GeV, Q 2 < 1 GeV 2 , and 0.3 < z < 0.9. The fusion production mechanism (curve 1), 
and the fragmentation production mechanism (curve 2) in the QMRK approach. 
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FIG. 5: pt spectrum of prompt J/ip in ep scattering with E p = 820 GeV, E e = 27.5 GeV, 60 
GeV < W < 240 GeV, Q 2 < 1 GeV 2 , and 0.3 < z < 0.9. The total contribution of the fusion 
mechanism (curve 1), the color-singlet part of the fusion mechanism (curve 2), the color-octet part 
of the fusion mechanism (curve 3), the total contribution of the fragmentation mechanism (curve 
4), the color-singlet part of the fragmentation mechanism (curve 5), the color-octet part of the 
fragmentation mechanism (curve 6). 
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10 2 




FIG. 6: z spectrum of prompt J/ip in ep scattering with E p = 820 GeV, E e = 27.5 GeV, 60 GeV 
< W < 240 GeV, Q 2 < 1 GeV 2 , and px > 3 GeV. The fusion production mechanism (curve 1), 
and the fragmentation production mechanism (curve 2) in the QMRK approach 
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FIG. 7: z spectrum of prompt J/ifi in ep scattering with E p = 820 GeV, E e = 27.5 GeV, 60 
GeV < W < 240 GeV, Q 2 < 1 GeV 2 , and p T > 3 GeV. The total contribution of the fusion 
mechanism (curve 1), the color-singlet part of the fusion mechanism (curve 2), the color-octet part 
of the fusion mechanism (curve 3), the total contribution of the fragmentation mechanism (curve 
4), the color-singlet part of the fragmentation mechanism (curve 5), the color-octet part of the 
fragmentation mechanism (curve 6). 
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